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Abstract 

We prove the existence of strong time-periodic solutions and their 
asymptotic stability with the total energy of the perturbations decaying 
to zero at an exponential decay rate as t — ► oo for a semilinear (non- 
linearly coupled) magnetoelastic system in bounded, simply connected, 
three-dimensional domain. The mathematical model includes a mechani- 
cal dissipation and a periodic forcing function of period T. In the second 
part of the paper, we consider a magnetoelastic system in the form of a 
semilinear initial boundary value problem in a bounded, simply-connected 
two-dimensional domain. We use LaSalle invariance principle to obtain 
results on the asymptotic behavior of solutions. This second result was 
obtained for the system under the action of only one dissipation (the nat- 
ural dissipation of the system). 
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1 Part I: Regularity and Stability for Solutions 
of a Magnetoelastic System in 3D 

1.1 Introduction 

We consider a system of partial differential equations that models the interaction 
between an elastic body and a constant magnetic field H acting on it. The 
elastic body is assumed to be conducting, non-ferromagnetic, homogeneous and 
isotropic. The body occupies a bounded, simply-connected domain £1 C K 3 with 
boundary of class C 2 . Let T be a fixed positive number and let u(x, t) denote the 
displacement field. Let f(x,t) denote an external force acting on the body. We 
describe the coupling between the elastic and the magnetic systems as follows. 
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The displacement in the body induces a magnetic field within the body given 
by H + h(x, t). The magnetic field exerts a Lorentz force on the elastic body to 
change the displacement field. The corresponding system of partial differential 
equations ([8] where p = 0) is: 

p M u" - //Ait - (A + /i) V div it + p(w') = /io(curlh) x (h + Hj +f, (1) 
h! + V\ curl curl h = curl u' x (h + H^j , (2) 
divh^O, (3) 

in Qt = X (0, T) (Qt = flx^ when we consider time-periodic solutions), u' 
is the partial derivative of u(x,t) with respect to t (time) and u" means («')'• 
Vector fields are represented in bold face. p M is the mass density per unity 
volume of the medium. A and \i are Lame's constants (/x > 0, 3A + 2/i > 0). 
v\ = l/(apo) where a > represents the conductivity of the material and /iq 
is a positive number representing the magnetic permeability. p(u') represents 
a dissipation that acts on the elastic body. 

The initial conditions associated with this system are 

u(x,0) = u (x), u'(x, 0) = Ui(x), h(x, 0) = ho(x), (4) 

while in the time-periodic setting, the periodicity conditions associated with 
this system are 

u(x,0) = u(x,T), u'(x,0) =u'(x,T), h(x, 0) = h(x, T). (5) 

The boundary conditions are 

m = 0, h-n = 0, (curlh) x n = on E T = dn x (0,T), (6) 

( in the time-periodic case: St = dQ x ^) where n = denotes the outer 
unit normal at x = (xi,X2,x^) G dQ. 

System (HJ)— (0) consists of a second order hyperbolic problem coupled with a 
parabolic equation similar to the Navier-Stokes equations through a semilinear 
coupling. 

The first investigation on the existence of solutions (with p = 0) for the 
magnctoclastic system in bounded domains was carried out by Botsenyuk [2], 
who proved the existence of weak solutions for the initial-boundary value prob- 
lem in £1 x (0, T) under the initial and boundary conditions ((4]) and ([6]). In 
[4], Botsenyuk proved the existence and uniqueness of a global strong solution 
under the hypothesis that the data is sufficiently small. 

Andreou and Dassios ([T]) studied the Cauchy problem in K 3 and showed, 
under some regular conditions on the initial data, that the solution of the system 
decays to zero, at a polynomial rate, as time goes to infinity. 

Menzala and Zuazua ([H]) proved the asymptotic stability for the system 
with linearized coupling and no mechanical dissipation (i.e., the only dissipation 
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acting was the natural dissipation of the system). The proof used LaSallc's 
invariance principle and did not provide a decay rate for the system. 

Rivera and Santos ([H]) proved for domains of special type that the energy 
associated to the linear three-dimensional magnetoelastic system decays poly- 
nomially to zero as time goes to infinity, provided the initial data is smooth 
enough. 

Charao, Oliveira and Menzala ([5 ) proved that the total energy of this 
system tends to zero as t — > oo when a nonlinear dissipation p(x, u') is effective 
on a small subregion of the domain. The rate of decay is given and depends 
on the behavior of p with respect to the second variable: algebraic decay for 
strong solutions and exponential decay if the behavior is close to linear for weak 
solutions. 

Mohebbi and Oliveira ([13]) proved the existence of weak time-periodic so- 
lutions for this magnetoelastic system under the effect of a nonlinear dissipation 
p(u ) like |i*|jgU with p £ [3, 4] (or the simplest case p = 0). The solutions have 
same period as the given time-periodic body force i(x,t). In that paper it is 
also argued that the existence of weak solutions holds for p > 2 (or p — 0) in 
the two-dimensional case. 

In this first part of the paper, we consider the regularity and stability of 
time-periodic solutions for the semilinear magnetoelastic system obtained by 
[T3] under the effect of a linear dissipation p(tt') = au' or a nonlinear dissipation 
satisfying certain restrictive condition. This results were motivated by our desire 
to show asymptotic stability of time-periodic solutions for the magnetoelastic 
system. 

1.2 Function spaces 

We consider a bounded, simply-connected domain 51 of M 3 with dil a 2- 
manifold of class C 2 . C§°(fl) and C£°(n) are the space of smooth functions with 
compact support in D, and the space of smooth functions with compact support 
in M. restricted to il. respectively. The usual Sobolev spaces are denoted by 
W s > p (n), W^ifl), s > 0, with the simplification H S (Q) := W s < 2 (n), flg(fi) := 
W°'\ti). The norm of H s is denoted by || || s and H 3 (fl) indicates the dual 
space of i?g(n). More details on some definitions that follow can be found in 

We recall that 

H(div,n) := {v E L 2 (Q) : divv e L 2 (Q)}, 

is a Hilbert space and C^°(fl) is dense in this space. The trace map j n defined 
on Cq°(Q) by j n (v) — v ■ n|gn can be extended to a linear continuous mapping 
7„ from H(div, VI) onto H~ 1 / 2 {dVt) which maps v into v ■ n\gn. Let 

V := {he C^(U) : div h = 0, h ■ n\ dn = 0} , 

Hl(£l) is the closure of V in iJ 1 (fi), which can be characterized as 

Hl(VL) = {he H^Q) :divh = 0, h ■ n\ dn = 0} . 
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Also, L 2 (Q) is the closure of V in L 2 (Q), which can be characterized as 

L 2 jn) = {he L 2 {VL) : divh = 0, h ■ n\ dn = 0} . 

The Hilbert space L 2 {Vt) is equipped with the usual L 2 -norm, | | 2 , and L 2 -inner 
product, ( , ) 2 ; while H^{Q) is also a Hilbert space with inner product 

((h,h)) :— (curl h, curl h) 2 

and norm || || induced from this inner product (which is equivalent to the usual 
ii^-norm). For s > 0, we define if*(0) as the closure of V in H s (n). H*(n) is 
a Hilbert space equipped with the inner product of ii s (f2), denoted by (( , )) s . 
As usual, H~ s (fl) indicates the dual of H*(Q,). 

Let H(cui\,n) := {v e L 2 (0) : cur\v e L 2 (fl)}, which is a Hilbert 
space with the norm ||w|| jf( CU ri,n) : ~ {M2 + I curlt^ 2 } 1 / 2 . C^°(fi) is dense 
in H (curl; VL). The trace -f T defined on C^°(S!) by j T (v) = v x n\gn can be 
extended continuously to a linear continuous mapping j T from if (curl, f2) into 
H- 1 / 2 (dn). We also define 

U 2 = {w e L 2 a n H 2 {Q) : curl™ € ff (curl, Q)}, 

r i */ 2 

equipped with norm defined by ||t«||(7 2 = ■{ || w lljj( cul .i q) + I curl curl to | 2 J- 

1.3 Theorem on the existence of periodic solutions 

For simplicity in notation, we consider the existence of periodic solutions of 
period T > for the semilinear coupled system rewritten in the following form: 

u" + C u + p(u') = curl hx (h + Hj + /, 

ti + v x Ch = curl u' x (h + Hj , (7) 

where 

C( ) = -n A( ) - (A + n) V div( ), V(C) = H 2 {Q) n ^(n), 
£( ) = curlcurl( ), V(C) = H 2 (n) n ff*(fi). 

We will find a weak T-periodic solution for this system under the following 
assumptions: 

(HO) ft C K 3 is a bounded, simply-connected domain with boundary 9f2 of class 
C 2 . 

Forpe [3,4], (f =(p + 2)/(p+l): 
(HI) / belongs to C([0, T],L q (fl)) and /(0) = /(T). 
(H2) The continuous function p : M 3 — > M 3 satisfies: 
There exists a positive constant K such that 

(p(z),z) e >k \z\ p + 2 , Vzem 3 , 
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There exist positive constants r p and K\ such that 

\p(z)\ E <KM V E X if|^>V 

(| \e and ( , )e denote the Euclidean norm and corresponding inner product.) 

Now, we state the theorem of existence of weak time-periodic solution proved 
by Mohebbi and Oliveira (2012) ([13]): 

Theorem 1. Let T > be the period of the function f. We assume (H0)- 
(H2). Then, there exist u G L°°(0, T; i^(ft) with u' G L°°(0, T; L 2 (Q)) n 
LP+ 2 {0,T;LP +2 (n)) and he L°°(0, T; L 2 (Q))nL 2 (0, T; if *(fi)) such that (u,h) 
satisfies the following system: 

- / (u',(p) 2 r]'ds + / a/ (u,ip)nds + / (p(u'), ip) 2 r]ds 
Jo Jo Jo 

= / B I (h,h)(<p)ri ds + / (/, ¥>) 2 ryds, 
Jo Jo 

y v e Hi (n) n l p+2 (n) , VijeD T , (8) 

- / {h,if>) 2 fj'ds+ / aii{h,ij))ds = / B H {u' ,h)(tf>)fjds, 
Jo Jo Jo 

VV G f^ /2 («), Vr) G D T , (9) 

w/iere 

D T := {cj G C°°(E) : w(s) = u(s + T), Vsel}, 
au{h,w) := v x ((h,w)), Wh,w G i^(fi), 
ai(u,v) :— fi (\7 u,\7 v) 2 + (X + /i) (div it, div r>) 2 , V«, f G ff (O), 

B u (w,b){ ):= J WX (&+#) • curl( )dcc 

Bi(h,b)( ) := y (curlfc.) X (& + #) • ( 

icitfi 5/j : L 2 (ft) x fl^(n) ->■ Ha 3/2 (fl) and Bi : H\{fl) x i?y 2 (f2) -> ff-^fi). 
Furthermore, u G C(0, T; i/ 1 / 2 (il)) and ft, G C(0, T; f7^ 1/4 (ft)). 

1.4 Strong Solution for the Initial Boundary- Value Prob- 
lem without Mechanical Dissipation 

We consider now the initial-boundary value problem (IBVP) described by the 
following scmilinear magnetoelastic system: 



u + Cu = curl h x ( h + H 

h + vi£h = curl 
div h = 0. 



u x (h + H 



(10) 
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1.5 A Priori Estimates 

We consider a priori estimates for solutions of the equations in the Faedo- 
Galerkin approach. For this we consider a basis {bjjjgN of eigenfunctions of 

ai(bj,ip) = yJa.Vy e i?o( fi ) 

and a basis {b^jfegN of eigenfunctions ([9J, P- 58) of 

au(b k ,ip) = A fe (b fe ,-0) 2 ,W> e J7*(fi). 

Given to € N, we define the finite-dimensional spaces S m = span &2j " " " ) ^m} 
and 5 m = span {6x, &2? • * • , &m}, so that 

m m 
^ ^ 6j , h in Cj bj . 

3=1 3=1 

We consider the weak formulation written in terms of {u m } and {h m } in the 
standard way. 

Let us define the energy of the system as 

2£{t) :=|u m (t)|2 + || Um (t)|| 2 + |h m (t)|2. 

It follows from the weak formulation with cp = u m and ip — h m , adding the 
resulting equations, the basic energy estimate 

§ + ^i||h m || 2 <|f| 2 |u m | 2 (11) 

This estimate implies that 

VW)< V^O) + ll/ll^(o,Ti£»(n)),VtG [0,T] 

using arguments similar to the ones that we will show for the second energy. 

It also follows from the weak formulation with ip — Cu m and i/> = Ch m , 
adding the resulting equations, the following estimate 

j t {||u m || 2 + |£u m | 2 + ||h„J 2 } + 2v l \lh m \ 2 < 2c, |£h m | 2 ||u m || 

+ 2||f||||u m || (12) 

where c, is a positive constant. 
Therefore, 

£i(t) + vi \Ch m \ 2 ds <£!{{)) + V2c, / £l l2 {s)\Ch m \ 2 ds 
Jo Jo 

+ V2 [ \\f(s)\\ E^sf^ds (13) 
Jo 



G 



where 

2S.it) := ||u m (f)|| 2 + \Cu m (t)\ 2 + l|h m (t)H 2 . (14) 
Brezis inequality implies that 

£l' 2 {t)<£l /2 {0) + V2c^ \Ch m \ 2 ds + V2 / ||f(»)||ds 

Jo Jo 

<c + ci f \Ch m \ 2 ds + c 2 [ \\f(s)\\ds (15) 
Jo Jo 



Returning to inequality (fT3|) . we get 

1/1 / \Ch m \ 2 ds<c 2 +c c 1 f \Ch m \ 2 ds + cj( f \£h m \ 2 ds 
Jo Jo \Jo 

+ 2c lC2 ||f(s)|| d s ^J Qf |£h m | 2 d^ 

+ c c 2 ^ ||f(fl)||dfl + ^n ||f(a)||<fej • (16) 
Let Z := J Q * |£h m | 2 ds. Then, we must have 

< c\Z 2 + jc ci + 2cic 2 ^jf ||f(«)|| dsj - i/ij Z 

+ c c 2 jf* ||f( S )||d S + ci(^l|f toll ^) + 4 (17) 

We recall also an auxiliary lemma due to Botsenyuk's ([!]). 

Lemma 2. Assume that 1,76 C°([0,T];R) are non-negative functions satis- 
fying the inequality 

x{t) < ~f(t)+ax 2 (t),Vt € [0,T], (18) 
where a is a positive constant. Suppose that the condition 

1 - 407(f) > 
holds for all t € [0,T] and £i(f), £ 2 (f) are roots of 

az 2 -z + 1 {t) =0, (19) 
6(f) < &(*)■ r/ien, j/ x(0) < Ci(0), it follows that x{t) < £i(f),Vf G [0,T]. 
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Using this Lemma, we assume that 

< c ci + 2cic 2 ( J ||f(s)|| ds 




4c?|c c 2 ^ ||f(s)||d« + c 
which simplifies to 

6<^fi(0) + 2i/ 1 V2c M f 1 1/2 (0) + 8c M ^i ( / ||f ||ds ] < ^. (20) 




Therefore J Q T |£h| 2 <is is bounded (by a constant that does not depend on 
m), which by ()15[) implies that there exists a positive constant Ce (which also 
does not depend on m) such that 

£i(t) < C E 

for t E [0,T]. These constants depend on the initial energy £i(0) and on the 
following norm of /: ||/||ii(o,r ; iji(n))- 

Remark 3. We conclude that if (|20l) is satisfied, there exist strong solutions 
for the three-dimensional initial boundary-value problem. This result improves 
on a previous result due to Botsenyuk (140 on the existence of strong solutions. 
Botsenyuk (TS) proved under his assumptions that there exists a unique strong 
solution. Under our assumptions the same result holds with the same proof. 

In the next theorem, the weak solution we are referring to are the ones 
proved by Botsenyuk ([2]). 

Theorem 4. Let T be a positive real number. Let f £ L 1 (0, T; ff 1 (51)). We 
assume that the following condition holds: 

6^(0) + 2is 1 V2c fi £ 1 1 /2 (0) + 8c^t ||/|Ui(o,r ; ifi(n)) < v\, (21) 

where 

2£ l {t) :H|u'(£)|| 2 + |/:u(t)| 2 + ||h(£)|| 2 , 

Then, the weak solutions (u, h) have the additional regularity: 
h £ L 2 (0, T; D{C)) n L°°(0, T; H^ty), 
u £ L 2 (0,T; D(C)) and 
u'gL»(0,Tiff{(fl)). 

1.6 Existence of T-periodic strong solutions 

We claim that this result of regularity holds for T-periodic solutions if the 
system has a linear mechanical dissipation (or nonlinear in the restricted sense 
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discussed later) , because we can show the existence of T-periodic solutions using 
the Poincare map. 

Using the Faedo-Galerkin framework, we take the L 2 inner product of each 
term in the first equation of the approximate problem (with p(u m ) = au m ) 

u' m + Cu m + p(u' m ) = curlh m x (h TO + H) + f 
by u' m , and the second equation, which is 

h' m + viCh m = curl u m ' x (h m + H) , 

by h m , we obtain the fundamental identity 

+ {p{u m ),u m ) 2 + ^i||h m || 2 = (f,u m ) 2 . (22) 



d£_ 



Formally, taking the L 2 inner product of each term in the first equation by 
eu m (with e > 0) and adding the result to the previous identity to obtain the 
following inequality 



d_ 
~dt 



\e + e(u m ,u m ) 2 + ^|u m | 2 } + ( 



'2 ~2~ 

12 , „ „ u ||2 < | f | .' 



rn 



" - £ >\ U m\ 

f| 2 |u m | 2 + C||h m || 2 



Let G £ := £ + e (u m ,u m ^ + ^|u m | 2 . Since 



-£ < G e < (2 + a)£ 



if < e < min(l, C n 2 ), we infer that 



(23) 



-G e + 2s£< |f| 2 |u m | 2 + C^ /2 |f| 2 ||u m || + C||h m || 2 \\ur, 
< C\{\ 2 V£ + C'V£\\h m \\ 2 
if < £ < min(l, C^ 2 , a/2,vi). This estimate implies that 



d 



at 2 + a 



^Gl{t) <C(|f| 2 + ||h m || 2 ) 



if < e < min(l,C 2 ,a/2,^i). 
Therefore, 



< C\ exp 



I e 2 +< 
Jo 



e^\f\ 2 (s)ds+ / e*J||h m || 2 (s)ds 



+ y/G E (0)exp 



2 + a 



(24) 



(25) 
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So, 



V£(t) < Ci ||f ||li ( 0,T;^(O)) + ^£ (0) (1 + ||f||zi(0,T;£ a (n))) 



+ ^llf 11^(0,^(0)) + (V2+^) V^(0)cxp f- 



■^-*) . (26) 
2 + a J 

where we have used the estimate (obtained from the uniform bound on \J~£): 

vi f ||h m (s)|| 2 d s < £(0) (1 + ||f \\ L i ( o,T;Vm) + H f \\h(o,T;mn)), V* e [0, T]. 
Jo 

We conclude that the S(u m (0), u m (0), h m (0)) = (u m (T), u m (T), h m (T)) maps 
every ball of radius R £ (i? cr , 1) and center at into the same ball, where 

C*i||f|| L i (0 ,T;L2 ( n)) + §f ||f|||i (0jT;L 2 (o)) 

-Kcr = — — 7 \ [4 I ) 

1 - x/2+Wp (-gfeTj - f (1 + ||/|U 1( 0,t^(O))) 

under the condition that i? cr £ (0,1). This map is also continuous. By the 
Brouwer fixed point theorem, it follows that there exists at least one solution 
(u m , u m , h m ) T-periodic. The remaining task of passing to the limit as m — > oo 
was described in detail in |13j and so we omit this part of the proof. The main 
conclusion that we obtain is the existence of a strong T-periodic solution when 
the dissipation is given by p(u ) = au , under the assumption that < R cr < 1. 

Notice that this result requires not only the expected restrictions on the sizes 
of i^i (large) and f (small), but also T not too small. 

Theorem 5. (Existence of strong T-periodic solutions) 

Iff£ L 1 (0,T;.H a (f2)) and the conditions ([2T|). < R cr < 1 hold, then there 
exist strong T-periodic solutions for the magnetoelastic system ([1]), ([2]), ([3]), ([5]) 
and ([6|). 

Remark 6. The drawback of this argument of existence depending on the reg- 
ularity estimate for the initial boundary-value problem and the Poincare map is 
that we have not been able to extend this argument for the full case considered 



1.7 Stability of Strong Time-periodic Solutions 

We proceed to investigate the stability of strong T-periodic solutions for the 
magnetoelastic system subjected now to a linear dissipation p(v ) = ay , For 
this type of dissipation we can use the result of the previous subsection since 
we proved the existence of T-periodic regular solutions using the initial value 
problem and the Poincare map. 
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Let (v, b) denote the perturbation pair of the strong solution (u, h). Then, 
these perturbations satisfy the following system 

v + £v + p(u + v ) — p(u ) = curlh x b + curlb x h + curlb x b 

b' +vxCh = curl[u' x b] + curl[v' x h] + curl[v' x b] (28) 
div b = 

under the initial conditions v(0) = v , v (0) = vi and b(0) = b . 

We will look for an a priori estimate for the solutions of this system. For- 
mally, taking the L 2 (S!)-inner product of the first equation by v , taking the 
L 2 (f2)-inncr product of the second equation by b and adding the results, we 
obtain 

~ {|v'|i + N| 2 + |b| 2 } + (p(u + v ) - P(u ), v ) 2 + ^||b|| 2 
= (curlh x b,v') 2 + (u x b,curlb) 2 (29) 

Let 2£ p (t) := |v | 2 + ||v|| 2 + |b| 2 denote the energy of the perturbation. Then, 
it follows that 



d£„ 

■ + \^p(u + v ) - p(u ), v + fi||b||* 

+ ||u|| ||b|| 3 / 2 |b| 2 /2 (30) 



^+ (p(u'+v)-p(u),v') 2 + z/i||b|| 2 < |£h| 2 ||b|||v | 2 



12 

Thus, using Young's inequality we infer that 



| + (p(u + v ) - p(u ),v ) 2 + |||b|| 2 < l|£h| 2 |v | 2 + 1 ||u || 4 |b| 2 



<-|£h| 2 |v| 2 + -C 2 |b| 2 (31) 



Assuming that V\ > 2\/2 CeCq , it follows that 



^ < -|£h| 2 |v'| 2 < -\Ch\ 2 2 E p (32) 
dt ~ vi 121 12 ~ vi 12 p V ; 



i.e., 



1 < —\Ch\ 2 iJTp (33) 



dt v\ 



which implies that 



V^< ex P [l l^h| 2 d S | < Chft^O), (34) 

since the term {/* |£h| 2 ds} is bounded by a constant as a consequence of the a 
priori estimates that produced the strong solution (u, h). This a priori estimate 
places the perturbation (v, b) in the usual energy space of the weak solutions. 
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Now, we take the L 2 (f2)-inner product of the first equation of (|28|l by ?yv, 
where 77 £ (0,a), to get 

»7^(v',v)2+7?(p(u' +v') -p(u'),v^) 2 + 77||v|| 2 -77|v'| 2 

<2»,||v|| ||b|| ||h||+T7||v|| ||b||||b|| Hl/2 

Adding each member of this equation to the corresponding members of equation 
(j3"Tj) . we obtain the estimate 

-^{£ p + r?(v',v) 2 }+7? (p(u + v') - p(u'),v) + (p(u + v') - p(u'),v' 
-^|v| 2 + r,||v|| 2 + ^||b|| 2 

< 1 |£h| 2 |v | 2 + c ||u'|| 4 |b|| + 2tj||v|| ||b|| ||h|| + t}\\v\\ ||b|| 3 / 2 |b| 2 /2 (35) 
^1 

We infer that 

j t {£ p + t}(v , v) 2 + ^|v| 2 } + (a - ry) |v | 2 + „||v|| a + ^-||b|| 2 
< - |£h| 2 |v | 2 + c ||u || 4 |b| 2 + 2r?||v|| ||b|| ||h|| + „||v|| ||b|| 3 / 2 |b| 2 /2 (36) 

Let 

G p ^ p + 77(v,v) 2 + ^|v| 2 . 
The next lemma establishes the relation between the function G p and the energy. 
Lemma 7. If < rj < min(l, (Cn)~ 2 ), then 

\£ P {t)<G p (t)<{2 + a)£ p {t). (37) 



Proof. 



l£n<G <£ v + l\v'\l + W 2r,2G s i ' m 



p i <-p T -] "|2 T ri2 2g2 1" " 2 _ l v l2 
< It/H + J|6|l + ( I + ^T=2 + ^ ) a '(^ «) < (2 + a)£ p . 



□ 

Then, 

^ + ( a ^)|v?+r / ||v|| 2 + ^||b|| 2 

< - |£h| 2 |v | 2 + c ||u || 4 |b| 2 + 2f,||v|| ||b|| ||h|| + „||v|| ||b|| 3 / 2 |b| 2 /2 (38) 
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So, for 1] e (0,a/2) sufficiently small, 

^ + f|v|^f||v!| 2 + |^l- C Ci||b|^ 

+ - 277 C £ } ||b|| 2 < „||v|| ||b|| 3 / 2 |b|f + 1 |£h| 2 |v'| 2 . (39) 

This inequality and the previous uniform bound on the energy of perturbations 
imply the following estimate 

^ + f |v'|| + |||v|| 2 + - c Clj |b| 2 

+ {^-2 v C E } ||b|| 2 < C hy /2£ p (0)r, £ p + A |£h| 2 £ p (40) 



Thus, 



^ + (c - C h ^2£ p {Q) ^ £ p < A |£h| 2 £ p (41) 



Now, using the previous lemma, 



dG p C -C hy /2£j$)r, 4 - 2 

+ — (2~+^) — Gp - ~! m2 ° p {t) [ ' 

Let 

_ C -C h y^(0)77 

Ci: ~ 

and 

C h := sup |£h|^dt,i > 
which exists by the previous estimates. Then, we infer that 

G p (t) < G p (0) exp (-dt + — 

This inequality implies, for 



< i] < min jl, ^2 • 

and 



2C 2 'C, V ^(0) 



^ > max | \/2 C B , 2c C% } 2C,- 1 / 2 , (43) 
the following main result 
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Theorem 8. (Unconditional Asymptotic Stability) 
If the conditions (|2ip . < R cr < 1 and (|43|) hold, then 

(i) there exist weak solutions for the initial-value problem (|28[) for the perturba- 
tion field (v, b) , 

(ii) the energy of the perturbations has an exponential decay to zero as t — > oo; 

£ P (t) < 2 (2 + a) £ p (0) exp (-C x t + . 

The same result holds for the case with nonlinear dissipation p(u ) (with K c 
replacing a ) under the assumption that there exists a positive constant K c such 
that 

{p(u + w) — p(u), z) E > K c (w, z) E , Vu, «i,z£l 3 



2 PART II: Asymptotic Stability for a Magne- 
toelastic System in 2D 

2.1 Introduction 

We consider now the asymptotic stability of solutions of an initial bound- 
ary value problem that models the interaction between a conducting nonferro- 
magnetic homogeneous isotropic elastic body that occupies at an initial time a 
bounded, simply-connected domain il C M 2 with smooth boundary. We suppose 
that a time-independent magnetic field (0,0, B ) with B : ft —> M. affects this 
region. Let u = (ui,u 2 ,0) denote the displacement field and (0, 0, h(xi, X2)) 
denote a disturbance of the magnetic field. Let T be a fixed positive number. 
We also denote by u the partial derivative of u(x, t) with respect to t (time) 

and u means ^ . Analogous notation is used for the scalar fields. 

The corresponding system of partial differental equations is the following: 

PM ~Q^ ~ /iAu ~ ( A + ^)Vdiv u + [i (B + h) Vh = f 2 (44) 
^-- Vl Ah + div ({B +h)^)= h (45) 

in Qt = O x (0,T). pm is the mass density per unity volume of the medium. A 
and p are Lame's constants (p > 0, A > Q). v\ = \/{apo) where a > represents 
the conductivity of the material and po is a positive number representing the 
magnetic permeability, f = (/i,f2) is a known external force. The initial and 
boundary conditions associated with the system are 

du 

u(x,0) = u (x), —(x 7 0) = u 1 (x), h(x, 0) = h (x) in il (46) 

and 

u = 0, Vb ■ n = on S T := dfl x (0, T), (47) 
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where n = (m, n 2 ) — n(x) denotes the outer unit normal at x = (xi,x 2 ) G dfl. 
The total energy associated with the above system is given by 



1 



2 



2 



£{t) = - ip M u r + M |Vu|i 2 + (A + ^Idiv u|£ 2 + |ft|£ a (48) 



L- 



where | and ( , )^2 denote the L 2 (f2)-norm and related inner product, 
respectively. 

If / = 0, then, formally, by taking the L 2 -inner product of each term in 
equation (|4"4"|) by u and each term in (|4"5j) by b and adding the resulting relations, 
we obtain the energy identity 

±£(t) = -v 1 \Vhf (49) 

which means that the model (T4"4"l) - f4"5)) has a natural dissipative term. This 
identity also provides the first a priori estimate ([3]) for each t G [0,T): 

\h(t)\h + f |Vft(t)|| a ds + |u (t)| 2 2 + ||u(t)|| 2 < Cl (50) 
Jo 

Botsenyuk ([3]) proved the existence and uniqueness of global solutions for 
the two-dimensional problem without assuming any restrictions on the initial 
data. The following estimate was obtained: 



ft|| 2 d S + ||ft|| 2 + ||u|| 2 + ||u|| 2 <c 2 , (51) 



where || || s denotes the norm of H s (£l), s G R, and < 7 < 1/2. 

We investigate the large time behavior of the solutions of this system without 
any additional damping mechanism (and no external force). 

2.2 Weak formulation 

Let T > 0, B n G iflfaTiH^fl) n L°°(Q)), h G i? 7 (^), u G # 1+7 (^), 
ui G fi G L 2 (0,T;L 2 (fl)) and f 2 G L 1 (0, T; iP(f2)) (0 < 7 < 1). We 

look for functions (u, u , ft) in 

L°°(0, T; H%(Sl)) x L°°(0, T; L 2 (ft)) x [L 2 (0, T; ff^fi) n L°°(0, T; L 2 (ft))] 

satisfying the following system: 

d f 
PMT(u'i) +a 2 (u,ip) +fi (B +h)Wh-ipdx = (f 2 ,ip), (52) 
* Jn 

W G Fq 1 ^) 

— (ft, y>) + ai(ft, ip) - (h, if) - j (B + ft) u' ■ V<^ = {fi,<p), (53) 
V^G ff 1 ^) 
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where 

ai(<p,rj) := vi (V^,V?7) + (93,77), V<p, 77 € if 1 (ft) and 
02 (v, w) := fi (Vv, Vw) + (A + div v, div w), Vv, w £ ifp (O). 
We will make use of the following existence-uniqueness result of Botsenyuk: 
Theorem 9. (Botsenyuk, 

Under the previous assumption on the initial data and forcing junction, then 
there exist unique functions (u, h) satisfying ^ , (|5|) , (|5^|) and (|53[) . Moreover, 
h £ i 2 (0,T;iJ 1+ T(^)), /i€C([0,T];^(f2)), u e C([0,T];ff o 7+1 (l])) and u' e 
C([0,T];iJT(fi)). 

2.3 Asymptotic behavior 

We shall assume that Bq is a constant and f = 0. We denote by L^{Vl) the se ^ 
of functions ip in L 2 (Vl) that satisfy / n </? dx = 0. We must have ho £ H 7 (il) n 
Lg(fi) and the solution h belongs to the spaces L 2 (0, T; n l|(n)) and 

C([0, T]; iJ 7 (51) nLg(il)), since (considering the three-dimensional problem) the 
identity dj(xkhj) — (djXk)hj, integration in f2, the boundary condition h-n = 
and the divergence theorem imply that /i must satisfy 

/ hdx = 0. (54) 
Jn 

We recall some definitions and results concerning dynamical systems. 

Definition 10. A dynamical system on a Banach space Z is a family 
{S(t)} t >o of mappings on Z such that 
(i) S{t)~E C{Z,Z),Vt>0, 
(it) S(0) = I, 

(Hi) S{t + s) = S(t)S(s)ys,t e R+, 

(iv) the function 1 1— > S(t)z belongs to C([0, 00), Z),\fz £ Z . 

An orbit (positive orbit) 7+ = 7 + (y) through ip in Z is defined to be "f + (ip) = 
U t >oS(t)ip. 

Definition 11. Let {S(t)} t >o is a dynamical system on a Banach space Z and 

V is a continuous scalar function defined on Z . The function V(<p) is defined 
by 

V(<p) = lim sup - [V(S(t)<p) - V(<p)] . 

V : Z — > M is a Lyapunov function on a set G C Z if V is continuous on the 
closure of G and V(ip) < for cp G G. 

Theorem 12. (JIUj) 

Let {S(t)} t >o is a dynamical system on a Banach space Z. LfV is a Lyapunov 
function on a set G C Z and an orbit 7 + (z) belongs to G and is in a compact 
set of Z , then S(t)(z) — > M as t — > 00, where M is the largest invariant set in 
S:={<p£G:V(<p)=Q}of{S(t)}. 
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Now, let X = Hl +1 {VL) x x [H~<{Q.) n Ll(n)] with < 7 < 1/2. Let 

F = Hl +1 (O) x fff'(O) x in'(n) n Lg(n)l with < 7' < 7 < 1/2. The total 
energy of the magnetoelastic system is defined by 

£(t) = \l [PM\u t (t)\ 2 + /i|Vu(t)| 2 + (A + M ) |div u(t)| 2 + fc 2 (i)] dz. 



We know that 



^ = -— / |Vfc| a dx<0, 



which implies 



£(*) + — / / \Vh\ 2 dxds = £(0), Vi>0. (55) 

°"Me JO JO 

Let 

f 7 (i):=||u,(i)|| 2 + ||u(t)|| 2 +7 + ||M<)|| 7 . (56) 

We also know ([3]) that there exists a positive constant C2, which depends on 
£(0) and also on ||ui|| 2 ,||uo|| 2 +7 , ||/io|| 7 , such that 

£y(*) + / ||/i|| 7+1 < C 2 , V<>0. (57) 



It follows from the well-posedness result of Botsenyuk ([3]) that {S(t)}t>o given 
by : X — »• Y, S(t)(uo, ui, /lo) = (u, u t , h) (t > Q), is a dynamical system. 

The continuous function V : X ->• R defined by V(z) = S(S(t)(z)), Vzel 
is a Lyapunov function for {S(t)} t >Q. 

It follows from (|55l) . (|57|) that for any initial data in X, the trajectory re- 
mains in a bounded set of X, for any t > 0. By compactness of the Sobolev 
embeddings, the trajectory remains in a compact set of Y. The hypothesis of 
Hale's theorem ([10]) are satisfied choosing Z = Y, G — X. Here, M is the 
largest (positively) invariant set in 



S = {(v,w,b)eX : [ \\7b\ 2 dx = Q}. 
Jn 



So, 

M= {(v,w,6) e X : (u(t),u'(t),h(t)) = S(t)(v,w,&) € X, 



/ |V/i(i)| 2 = 0, Vt > 0}, 
Jn 



i.e., 



M = {(v, w, b) e X : (u(i), u (i), fc(t)) = S(t)(v, w, 6) € X, 
h{x, s) = h(s) a.e. in Q, Vi > 0}. 
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Choosing ip = 1 in (I5U1) . we conclude that ft- = constant. It follows from ([54")) , 
that ft = 0. 

Since we are assuming that Bq is constant, it follows that M is given by 
M = {(v,w,6) € X : (u(t),u (*)»M*)) = S(t)(v,w,6) € X, 
ft(x, s) = 0, div u = 0, pmvl + Cu — in il x [0, oo)}. 
If, in addition, f2 has the following property V 

If £ € Fo(fi)is such that 

(i)-A£ = 7 2 £, (w)div£ = in ft, 

(m)£ = on dfl (58) 

for some 76R, 

imply that £ = 0, 

then, M = {(0,0,0)}. 

We have proved the following result on the asymptotic behavior of the system 
under the action of only the natural dissipation that comes from the equation 
for the magnetic field: 

Theorem 13. Let Q be a bounded, simply connected domain in R 2 with bound- 
ary of class C 2 . Let < 71 < 72 < 1/2. Let {S(t)} t >o be the dynamical system 
defined on X 7l = # 1+71 (fi) x (O) x [H^ (fi) n L 2 (f})] &2/ £(i)(u , Ui, ft. ) = 
(u(f), Ut(f), /i(t)) (i > Oj. We assume that Bq is constant. Then, 
(i) If has the property V , then 

u^0 (in F 1+7l (O)), ft->0(miJ 7l (^)n^(r2)), as f 00. 

(n,) // does not satisfy property V (e.g., is a disk), then for any initial 
data in X l2 , (u, u t , ft) (t) aproaches the set M in the norm of the space X lt as 
t — > 00, i.e., 

u^v(mff 1+7l (O)), ft^0(mi? 7l (fi)ni^(^)), as t -> 00, 

where v satisfies div v = 0, v + Cw = 0, where Cu := — /iAu — (A + div u. 

Remark 14. Related to the property V are the following facts: 

(i) Weinacht (1151) proved that if fi is a simply connected two-dimensional 
bounded region with boundary of class C 2 and the eigenvalue problem 

Czi = Xi zi in O, (59) 
div z; = in f2, (60) 
z;|on=0 (61) 

has classical non-trivial solutions for each I — 1, 2, • • • , then f2 is a disk. 

(ii) In Dafermos (fSjl), it is observed that for a disk D = {x G R 2 : \x\e < 1}, 
there are non-zero eigenf unctions v m — (w™,w™ = (0, Ji(Cm r ))> rn&N, where 
Ji is the Bessel functions of the first kind of order 1 and £ m is the m-th non- 
negative root of J\ . 
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